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Abstract 

Expressions for the thermodynamic potential of a Dirac fermion gas 
are represented at finite temperature with the chemical potential in an 
ultrastatic space R d x S N . The high- and low- temperature expansions 
for the thermodynamic potential are obtained and, in particular, strongly 
degenerate fermi gas is investigated. For the Candelas-Weinberg model, 
sufficiently high "charge" density prevents the compactification of the ex- 
tra space. 

1 Introduction 

These days many efforts have been made for pursuits of "the unification through 
higher dimensions." It is highly remarkable that supergravity theories have 
simple or unique structures in more than four dimensions. [1] Recent success in 
superstring theories suggests that higher dimensions are required not only for 
simplicity, but also the consistency of the theories. [2] Anyway, as far as we regard 
the extra dimensions as physical settings, we obviously need to comprehend the 
phenomena called "compactifications" , in order to find out our four-dimensional 
world. 

Recently, there appear many works investigating the behavior of the time de- 
pendence of the scale factor of extra spaces, that is, the "Kaluza-Klein cosmology." [ 
Many people intend to explain the large amount of entropy in our universe in the 
same scenario. The scenario is referred to with a slogan "entropy comes from 
extra dimensions" [4] for this attempt. Therefore it is necessary to calculate the 
thermodynamic quantities in higher dimensional as well as curved spaces. 

On the other hand, Actor derived the thermodynamic potential in arbitrary 
dimensions. [5] Particularly, we are interested in the case with non- vanishing 
chemical potentials. It is well known that symmetry restorations depend on 
finite density effects as well as finite temperature effects in four-dimensional 
universe. [6] Thus, one can expect that compactifications of the extra spaces are 
similarly influenced by finite density and temperature. [7, 8] 
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Since we want to know, at least, about the thermodynamic quantities before 
and after the compactifications, we have to obtain the thermodynamic potentials 
in curved-spaces. In the present paper, we derive the thermodynamic potential 
of a Dirac fermion field in curved spaces such as T x R d x S N . 

The present paper is arranged as follows. In §2, we briefly review the ther- 
modynamics and the effective potential in field theories with imaginary time 
formalism. [9] We treat only with a Dirac field throughout this paper. In §3, 
we show the high temperature expansion of the thermodynamic potentials in 
the space whose background metrics are that of T x R d x S N . The low tem- 
perature case is studied in §4. The strongjy degenerate femi gas in higher- 
dimensional and/or curved spaces arc studied. In §5, we discuss the "Kaluza- 
Klcin thermodynamics," [10] especially in the case of non-vanishing chemical 
potentials. The last section is devoted to a summary and discussion. 



2 The thermodynamic potential and the path 
integral 

We introduce the thermodynamic potential il for an ensemble at temperature 
/3 _1 and chemical potential /x as follows: 

Z G = exp(-^O) = Tr cxp{-/3(iJ - fiN)} . (1) 

Here H is the Hamiltonian of the system and N is the particle number operator. 
Zq is the so-called grand partition function. Equation (1) reads for fermions, 
after the mode expansion as 

Z G = exp(-/3f2) 

= JJ[exp(- J 8w fc )(l + cxp(-/3(w fe - /i)))(l + cxp(-/3(w fe + /x)))] , (2) 
k 

where the frequency of modes is given by u>k = (k 2 +M 2 ) 1 / 2 . Here the zero-point 
oscillation is included. 

Now the entropy S, the pressure P and the particle number N for the system 
under consideration are given by the partial derivatives of f2(/3, V, /x), 

r-fv 

The "particle number" introduced here represents "the particle number mi- 
nus antiparticle number," which may be called "the charge asymmetry" . 

In the path integral language, the thermodynamic potential for a Dirac field 
with mass M is related to the one-loop effective potential, [5] 

lnDct[#>+ M] = Trln[(w„ + i/i) 2 + ui 2 k ] with u? k = k 2 + M 2 , (4) 

where Tr means the integration and summation of all physical modes and degrees 
of freedom (i.e., including a trace of the Dirac matrix). Here, uj n = {2n/f3}(n + 
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1/2), (n is an integer) is introduced in the imaginary time formalism for field 
theories at finite temperature. [9] The chemical potential in (4) can be regarded 
as the zeroth component of (imaginary) gauge field which is coupled with the 
charge density of the field. [5, 11] 

We can find the relation between one-loop effective potential given by (4) 
and the "quantum mechanical" [8] expression of the thermodynamic potential 
(see (2)) through the following identities (for fermions):[12] 



{n + 1/2) + in) +y 



= (5y + ln(l + exp/3(^ - y)) + ln(l + exp -p(p + y)) . 



(5) 



It is well known that the quantum vacuum energy at zero temperature 
comes from the first term on the right-hand side of (5). Apart from this zero- 
point energy, we derive the thermodynamic potential for a Dirac field in flat 
d-dimensional space using (5) (and expansions of the logarithms): [5] 



(t 



(2BM\ (d+1) ' 2 
x2(^j %i)/2Pn) 



(6) 



where K v {x) denotes the modified Bessel function, V d is the volume of d- 
dimensional space and trl = 2^ d+1 ^ 2 V 

We know another way of evaluating the one- loop effective potentials. In 
general, the one-loop quantum corrections contain divergences which need to be 
regularized. Using the zeta function regularization,[13] we obtain the following 
expressions for effective potentials: 



/^ = C'(0)+ln(2^4)C(0), 



where 



as) 



trl V, 




+ k + M 



(7) 



where fi R is a parameter which has the dimension of mass and comes from ad- 
justing the scale of the measure of the path integral. [13] Very recently, Allen 
showed that p? R is only involved in the vacuum energy part of the thermody- 
namic potential which is independent of /3. [14] Therefore, as far as we discard 
the vacuum energy, we do not need to worry about the regularization-scale [i R 
and the thermodynamic potential of the system can be derived from the first 
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term of (7). To see this, we use the following identity (see the Appendix): 

1 + 2 ^ (-1)" cosh(n/3 M ) exp ( - ^- J 

n=l ^ ' 



P ,-1/2 



(47T)V2 



(8) 



Then we can divide ((s) into two parts: 



as) = 

Co 00 = 



trl g V d 
2 (4f)V2r(s) 
trl /3 V d 



2 (47r)V2r(s) 
xcxp|-t(k 2 + M 2 )-^| . 



r d " s " /2 /(^ cxp{ ^ (k2+M2)} ' 

l dtts ~ 3/2 J ( 2 ^ 2 E(- 1 )" cosh (^) 



(9) 



Co (0) contributes to the thermodynamic potential as the temperature-independent 
vacuum energy. Thus, we only consider the contribution from C/j(0)- I n f &c t, 
r(s)C / g(s) does not diverge when s — > and T(s) ~ s _1 , so we can easily find: 



^(0) 
trl V d 

~Y (4tt) 



1/2 



-3/2 



(27r) d 



x2^(-l)"cosh(n/3 A i)expi -t (k 2 + M 2 ) - 

■n — 1 ^ 



/3 2 n 2 



71=1 

tn y d 



-(d+l)/2-l 



2 (47r)( d +!)/ 2 7 

oo , 

x2^(-l)"cosh(n/3^) exp j -iM s 

n=0 



f3 2 n 2 
At 



trl 



(47r) (/ +1) /2 ^ (d+1) E(- 1 )" cosh (^) 



x2(^^-j K {d+1)/2 (/3Mn) 



(10) 



This coincides with (6). 



3 The high-temperature expansion 

In this section, let us derive the thermodynamic potential in curved space. We 
consider the static backgrouud metric of R d x S N as a simple but sufficiently 
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general case. In order to treat this case, we generalize (10) as follows: 
tr 1 V d 



n = 



2 (47r)( d + 



— \ dtt-V+^-^^i-ircosHnM 

J° n=l 



00 ( R% 2 



(11) 



where 



r(jv + i) 



zir(AT) 

and trl = 2l D / 2 l with D = d + N + 1. a is the radius of the hypersphere S . 

Now, we consider high-temperature case such as j3/a <C 1. In this case, we 
can use the following formula for infinite sum related with theta function to 
expand (11) with respect to (3/a: 

2 "1 



^ di exp • 



1=0 



, N 



1 T(N/2) 



.N/2 



1 - — N(N- l)x + 0(x 2 ) 



(12) 



2 T(N) 

This identity is explained in terms of theta function. (See the Appendix.) 
Making use of Eq. (12), we can expand Eq. (11) as 

V d V N 



Q = (tr 1) 



(4 7r )(d+Af+l)/2 

(2pM\ [d+N+1) ' 2 



p-(d+N+i) cosh(n/3^) 



n=l 



K {d+N+1)/2 (f3Mn) 



\ n 

N(N-l)f3 2 /2(3M\ (d+N - 1)/2 
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a 2 \ n J 



K (d+N _ 1)/2 (/3Mn) + • 



as) 



where V N = 2TT ( - N+1 ^ 2 a N /T((N + l)/2) is the volume of S N . This first term on 
r.h.s. of Eq. (13) is the thermodynamic potential in flat (d + iV)-dimensional 
space, while the second and further terms correspond to the deviation from the 
flat-space case. 

When fi = 0, one can and the same expression as Dowker's.[15] For even D(= 
1 + d + N), the summation over n becomes a polynomial (i.e., finite terms). [5] 
We show here the expression in the simple massless case for later convenience: 

n = (tri) ,, J^ 1W9 /r( d+JV+1 > 



(4 7 r)( d+Ar + 1 )/ 2 



2 d+AT+lp 



d + N + l 



£(-i)" 



( cosh(n/3/x) 



n=l 



JV-1t 



d + iV- 1^ N(N-l)/3 
12~ 



2 00 



^ cosh(n^i) 

> „d+N-l 



n=l 



Id) 
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4 The low temperature expansion and the strongly 
degenerate Fermi gas 

Actor gave the low temperature expansion of the thermodynamic potential. [5] 
However, the low-temperature approximation becomes meaningless when M 2 < 
y? . When M 2 < /i 2 , it is well known that Femi gas degenerates strongly at 
low temperature. Fortunately, we know the formula for the polylogarithmic 
functions Lijy(x) — X^^Li x n /n N :[16] 



LK-2T 1 ) = (-1) 

N-2 



N-l 



Liiv(-2/) 



E 

r=0 
til 

Nl 



(-iy 



r\ 



(lnyni-(-l) A '- 1 -n(l-2 



-,r-N+l 



)C(N - r) 



N 



-(In y) 



N 



(15) 



and this gives another method for expanding f2 at low temperature. Hereafter, 
we consider M = case, for simplicity. The generalization to the massive case 
is straightforward. 

In (i-dimcnsional flat space (which means dimension of space-time is 1 + d) , 
the thermodynamic potential for a massless Dirac field is given by (cf. (14)) 



n = (trl) 



V d 



trl 



(4 7r )(d+l)/2 

Vd 



^-(d+l)2<i+lp 



(d+l) 2 d+lp 



, cosh(n(3n) 



r d+i 



n=l 



( 47r )(d+l)/2^ ~ ' \ 2 

Applying (15) to (16), we can reexpress it as 
trl V d 



[Li d+1 (-e^)+Li d+1 (-e-^) 



n = 



2 (4 7 r)( d + 1 )/ 2 

(d + iy. 



[(d+l)/2] 



l + (d + l)\ Yl 2 ^ 2h ) 



-2b 



6=1 



{d + 1- 26)! 



(1 



•,l-2b\ 



+{l + (-l) d }^d+i(-e-^)} 



(17) 



It is apparent that this expression enables us to approximate itself at low tem- 
perature by evaluating the summation and the polylogarithmic functions ap- 
propriately. 

Furthermore, if d + 1 is even, we get the exact form of the thermodynamic 
potential which is expressed as polynomials: 

,d+l 



trl V d 



2 (4^)<V 2 r(^) d+1 

(d+l)/2 

1 + (<*+!)! £ 



-26 



6=1 



(d + 1 -26)! 



2(1 



,1-26 



)C(26) 



G 



trl Y± a-(rf+l)od+l r 

2 (4 7 r)(rf+ 1 )/ 2P 



E iir 2 ^^-^-^- d )+ ((i+1)! 

c=0 



(18) 



This is exactly the same expression as the one obtained by the high-temperature 
expansion. [5] 

Let us turn our attention to the space R d x S N again. Integration over t in 
(11) yields when M = 0, (cf. (10)) 



V 

n = (tn) ( ; +1)/2 r (rf+1) E(^ 1 )" cosh (^) 

^ ' n— 1 

°° /2/?lj,\ (d+1)/2 

x E( 2 ^) 2 Hr ^( d+ i)/2(^n). (19) 



Using the integral representation 



K v {z) = £pW I™ e-^-Xy-^d^ (20) 

r(i/ + 5) A 



the sum over n can be performed to give 



(21) 



Now we can analyze the low-temperature case by using the method which can 
be found in many textbooks. [17] 

First of all, we examine the zero-temperature (/? — > 00) case. Note, in the 
limit (3 — > 00, 

'(-*>• ^ 

where #(2) is the step function. The expression (21) reduces to the following 
form by means of (22) at zero temperature: 

« = - trl (47 ^/ 2r( L ) E - > ( 23 ) 

where Z TO is the largest integer satisfying < fi and if oj > then fi = 0. 

As a check on the efficiency of this approach, let us consider the limit a — > 00. 
In this situation, the sum over I is reduced to the integration by the substitu- 
tions: 

u^z, J2 d ^YW } J o dzz N -\ (24) 
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Consequently, we can show: 
(4' 

v d v N 



trl Vd - — f ' dzz N+d dx (x 



1) 



d/2 



-trl 



( 47r )(d+iV)/2 r (d±2) r ( ^) d + TV + 



trl V^Vat 



2 (4^)( d + Ar )/ 2 r (^^±f±^) rf + iV+l 



(25) 



As is expected, this result corresponds to the one obtained from (17) in flat 
(d + 7V)-dimensional space at zero temperature. 

Particularly, we are interested in the case d = 3. In this case, one finds: 



= -4 



Vz 
24tt 2 



5> 

1=0 



Mm 2 - "?) 1/2 U 2 - H 



2 |Wi W 



1/2^ 



(26) 



From this expression, we realize that (26) is merely the sum of the thermo- 
dynamic potential in four-dimensional space-time at zero temperature [18] over 
the discrete mass levels which arise from compactification of N dimensions. The 
higher-modes than the chemical potential are frozen out. This interpretation is 
easily understood by drawing the figure which is analogous to the one in the 
paper by Barr and Brown, [4] but where we replace T with \i. Of course, in our 
case, the extra space S N does not admit zero modes, therefore the situation is 
slightly different. If we consider the case that the extra-space has zero modes, 
we will obtain the exact expression for Q, in the flat d-dimensiollal space when 
fi is smaller than the first non-zero massive mode. 

Before considering finite but still low temperature effect we show another 
example, the thermodynamic potential in S N (i.e., d = 0), 



SI = -(trl) 2*0* 



1=0 



-(trl) 



(l m + N)l (\ia l m + (N + l)/2\ 1 



UiN-iy. \n 



N+l 



(27) 



where 



N 



By using relationship (3), we derive the particle number 



(28) 
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We can see the effect of the discreteness of the Kaluza-Klein mass level from (28). 
One can also find self-consistent solutions of Einstein equations after balancing 
between the stress tensor of the degenerate Fermi gas with the Casimir effect 
and the cosmological constant as the case at finite temperature considered by 
Dowker et al.[19] The cosmological constant in this case plays an analogous role 
to the bag constant in MIT bag model. [18] Although the above example may 
be an interesting exercise, we do not study this situation in this paper. 

Now, let us return to treating low-temperature effect applying the textbook- 
formula [17] to (21), we can get the following low-temperature expansion for 
fi: 



x M {(d - 1) M 2 - 3cjf}(^ W f) d / 2 - 3 + ...]. (29) 

For the manifold which has zero-modes in contrast with S N of our case, one 
can find that the expression for the thermodynamic potential agrees with that 
in flat d-dimensional space derived by (17) when < fi < (the mass of the lowest 
massive mode). 



5 Kaluza-Klein thermodynamics and instability 
at finite density 

In this section, We first examine the themodynamics in the space R d x S N . 
Similar case with vanishing chemical potential is extensively investigated by 
Tosa,[10] so we will concentrate our attention to the case at finite density. The 
remainder of this section is devoted to the investigation of instabilities of com- 
pactifications induced by quantum effects [20] at finite density. 



Kaluza-Klein thermodynamics at finite density 

We consider a massless Dirac field in the space R d x S N . under such circum- 
stances, we assume the thermodmamic potential f2 must be of the form 

0=I(f) d / ( ,a^/a).I(f) d / ( ,,,), (30) 

where R is the scale factor of the flat <i-dimensional space and a is the radius of 
the iV-dimensional hypersphere S N . 

The chemical potential is introduced for a fermion field in (d+7V)-dimcnsional 
space. In the Kaluza-Klein sense, we find indefinitely many particles in d- 
dimcnsional space, however, the particle number defined in this paper is to be 
conserved, because this may be regarded as the charge asymmetry in the system. 
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Now let us derive several thermodynamic quantities using Eq. (30). We 
obtain 





(31) 



where V d , Vn, P, Q, S, AT and E are, the volume of R d , the volume of S N , the 
pressue for R d , the pressure for S N , the entropy, the particle number and the 
internal energy, respectively. We notice immediately, 



This relation has often been mentioned, and is equivalent to state the thermal 
contribution to the stress tensor is traceless. 

On the other hand, the conservation of the stress tensor gives 



It is natural to assume the entropy and the particle number are conserved sep- 
arately. This point will be discussed again in the next section. 

Instability of compactification at finite density 

In order to gain stable compactifications of extra spaces, various mechanisms 
are proposed. Candelas and Weinberg [20] considered quantum effects of matter 
fields in the case of the compact spaces, S N . Many other authors showed the ex- 
tended versions of their model and the presence of various sorts of instabilities. [21, 
22] Recently, Accetta and Kolb exhibited the finite temperature instability for 
compactifications, and the critical temperature for the instability. [23] Here, we 
show the finite density instability, which can be discussed almost parallel to the 
work by Accetta and Kolb. 

At high density (/3p >• 1), the themodynamic potential for Dirac fields can 
be written in a form 



-p + d- P + N ■ Q = 0, here p = 



E 



(32) 



V d V N 



dE + PV N dV d + QV d dV N = TdS + fidM = . 



(33) 



{C N - C(pa) 



d+N+l 



) 



(34) 



a' 



d+i 
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Note here fi contains the Casimir stress energy which is determined at one-loop 
level when the space-time dimensionality is odd. The pressure for the internal 
space, Q, and the particle number N are computed from (34): 

Q ~ V N N a 8a 

] " ' [(d + l)C N + N-C^a) d+N+1 ], 



V N a d +! 

on 

Y± C {d + N +\){,ia) d+N 



(35) 



We recognize the similarity to the high-temperature case. [22] The only difference 
is the exchange of /j, and T. Therefore we can immediately obtain the critical 
value for jj,: 



f^crit 



1 

a 



C N 2{d-\) (d+1 



C 



N 



N 



l/(N+d+l) 



(36) 



where ao is the static radius of S N . Using this, we conclude that there are no 
stable compactifications when Af/Vd > (Af /Vd) C rit, where 



C 



d + N+l 



crit a 

For example, we take d = 3 and N = 7. In this case, 



d 



C N 2(d- 1) (d+1 



C N 



N 



+ 1 



(d+N)/(d+N+l) 



(37) 



C 



5.958744 x 1CT 5 x / , 

trl 1 2 7 r( jv + 1 )/ 2 /r (*±1) 

2 (4 7 r)(d+N)/2 r (d+N+2^ d + N + 1) 

3.139894 x 10" 7 x /, 



where / denotes the number of Dirac fields. Then we obtain 



'— | = 4 x 3 - 35 x 111 

•it a o 



(38) 



(39) 



If we consider the initial size of 3-space is the same order as the one of the 
compact space, Eq. (39) reads the severe constraint to the charge asymmetry. 

One can discuss more details (such as including semiclassical instability [22]), 
by repeating similar analyses of Accetta and Kolb.[23] 



6 Summary and discussion 

We gave the expressions for the thermodynamic potential of a Dirac fermion 
gas with the chemical potential in the space R d x S N . We mainly pay our 
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attention to the low-temperature case which leads to the strolls degeneracy 
of fermi gas. We also gave some examples for applications to a few aspects of 
Kaluza-Klein theories. We omitted the detailed discussion about the application 
to the Kaluza-Klein cosmology, which thus will be discussed in another occasion. 

The treatment of finite density implicitly assumes the existence of conserved 
U(l) charge. Superstring theories imply the presence of the gauge held as "pri- 
mary held," so we may expect some conserved charges. The chemical poten- 
tials accompanied with non-abelian charges are also introduced by Haber and 
Weldon.[ll] The generalization to the high-dimensional case will enable us to 
investigate the finite density effect on the breakdown of primary gauge sym- 
metries as well as compactihcations. It may also be interesting to consider the 
themal property of the gauge held itself [24] in higher dimensions. 

On the other hand, it is also interesting to consider the chemical potentials 
associated to each "Kaluza-Klein charge," which are induced after the compacti- 
acation. This analysis will become necessary when one studies the property of 
the pyrgons [25] at finite temperature and density. We leave this for future 
publications. 

We need to consider the thermodynamic quantities for boson fields for further 
study. There is the technical subtlety to calculate the themodynamic potentials 
in relation to the bosonic field with the chemical potential. [26] But there are also 
attractive phenomena such as the Bosc-Einstein condensation. We will report 
on these problems generalized to higher-dimensional case elsewhere. [28] 
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Appendix 

We explain the identity appeared in §§2 and 3 in terms of theta functions. First 
we observe: 



$3(u|t) = exp(i7TTn 2 + in2nv) , 



(40) 



and the well known relation: 



#3 




= (— it) 1 ! 2 exp(iirv 2 /t)$3(v\t) 



(41) 



Using this, (8) can be derived as follows: 




) 



12 



U 2 



It**) tf 3 (21(^+- 



L 

2tt 



47Tt 



47rf 



2 n 1/2 oo 



S exp - 



4i 



(4tt) 1 /2 



-1/2 



exp ^ ^z/i - 

pV 



l + 25>D"exp(-£JM cosh(Ai ^ n) 

n— 1 ^ ' 



(42) 



Next, we show the derivation of (3.2). A similar expansion but for bosons is 
given by Yoshimura.[27] We consider the case N is odd. 



S F = '^2diexp{-(l + N/2) 2 x}, where di 



1=0 



T(l + N) 
l\T(N) 



(43) 



We write the degeneracy di as 



di = 



T(N) 



V -~2 



V+ 2 



V ~2 



l + v + 



V ( 1 

= C » m ( 1 + v + 2 

m=0 ^ 



TV - 1 



C„m arc independent of N, and the first two terms are given by 



C vv — 



T{N) ' 



1 1 

24T(/V) 



N(N — 1)(N — 2) 



(44) 



(45) 



(46) 



Since the sum over to vanishes when I = — 1, — 2, • • • , — v, we carry out the 
summation over I from I = — v to oo. Then, substituting (44) into (43), we 
obtain 



Sf — ^ C vm {— 1) 

m=0 



dx m ^ 6XP 



1 



m=0 

Here, we have used the identity derived from (41) 

2 



n— 1 v 



.(47) 



H exp 



i = — oo 



1 + 



1 



x 



exp 



(-1) 



IX 

2^ 



ix\ 

7T / 
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- ^f(-i)-«p(-^)- w 

n— — oo x 7 

From (47) and (44, 45, 46), we find the asymptotic form of 5p, 

Sf m^*- N/2 (l - ^(AT - l)x + 0(* 2 )) . (49) 
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